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SUMMARY 



A comparison is made of the results of calculations 
of the critical flutter speeds of nine uniform rectangular 
wings without ailerons by two methods , as follows: 

1* The method in which the effect of the wing flutter 
shape is introduced indirectly by choosing aean values for 
t h e w 1 n g para ,m e t e r s 

2. A method that introduces an additional decree of 
freedom in the wing— "bending flutter node 

These examples show that, for the case of normal 
wings and even for wings of higher aspect ratio, the dif- 
ference in the values obtained with these two methods is 
inconsequential and the effect of including a second de— 
gree of freedom in the bending mode is usually small. 



INTRODUCTION 



The effect of the shape of a fluttering airfoil is 
int r educed into the equat ions of mot ion and a compar is on 
is ma.de of the predicted critical flutter speeds obtained 
by the several methods that have 0 eon developed f or s Giv- 
ing the flutter problem. 

Lagrange T s equations, applied to the flutter problem, 
pr ov i de a means of s Giving f or the cr it i cal f lut t er speed 
and frequency and, in principle, for the wing flutter shape 
In the present paper, flutter speeds predicted by the use 
of this method for uniform wings of rectangular plan form 
are compared with the results obtained by the standard two- 
dimensional method (reference l) in which the effect of 
wing flutter shape is only indirectly introduced. Soth 
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values are compared wi|fc«iia*4l4flp» speeds obtained from 
experiments made in the 8-foot high-sp^d tunnel at LMAL . 
This comparison is used to judge the importance of intro- 
ducing considerations of wing flutter shape into the 
flutter problem. 

A solution of the flutter problem was given by 
Th^odorsen (reference l) and its application and implica- 
tions were amplified by Theodorsen and G-arrick in a sub- 
sequent paper (ref ere nee 3 ) . The f or mal s olut i on of th* 
flutter problem was carried out on a uniform infinite wing 
about which the air flow is two-dimensional. The struc- 
tural stiffness of the wing in bending and torsion is 
represented by equivalent springs. The aileron is assumed 
to extend the ent ire length of t he wing end is supported 
at its hinge by a spring r uprose nt ing the elastic pr oper- 
ties of the aileron-control system. The wing has three 
distinct types of vibratory motion: bending motion per- 
pendicular to 'the plane of the wing, torsion motion about 
the static torsion axis, and aileron motion about the 
aileron hinge, Flutter is a sustained oscillation involv- 
ing two or more coupled vibration modes. Four types of 
flutter are thus recognized: (l) be nd ing-t or s i on ; 
(2) bonding— aileron; (3) torsion-aileron; and (4) bending- 
torsion-ail^r on. The primary parameters determining the 
critical flutter speed are found t o be : the wing chord, 
th^ mas s of the wing per unit span length , the chordwise 
position of the static torsion axis and center of gravity, 
the moment of inertia of the wing about the static torsion 
axis, the air density, the torsion and binding stiffness, 
th* aileron cent of —gravity position with r^sp^ct to the 
aileron hinge, the aileron hin^e position on the wing chord, 
the aileron mass per unit length, anc 1 the stiffness of the 
ailer on-c antrcl system. The mas s and s t if f ne s s parameters 
are combined in the equat ions of mot i on of the wins to give 
nondinens ional coefficients that define t he bending , t or- 
sion, end aileron natural vibration frequencies. 

APPLICATION OF THE TV0-DIM3NS ICiTAL FLUTTEH THEORY 

TO FINITE tflFGS 



When the two-dimensional flutter theory is applied to 
a wing that is a finite continuous structure with section 
paramo t er s var iable al ong the span , several que stions ere 
raised: 



1. What are the air forces assignable to a wing 

vibrating with amplitude variable along the 
span ? 

2. Which section parameters of the wing are to be 

used in solving for th^ critical flutter 

s p ■ ?. e d ? 

3. What values' are to V-> given to the frequency 

terms occurring in the theory? 

The present practice in applying the two—dimensional 
flutter theory is to .pick a position along the wing span 
that is considered -representative of the averrga proper- 
ties of th* wing with regard to section parameters and 
wing flutter motion* The three-quarter semispan section 
is normally considered to he this characteristic wing 
section. The lowest natural vibration frequencies in 
wing bending , torsion, and aileron oscillation are 
assigned. The solution for the critical flutter speed 
is then assumed to be that of a uniform infinite wing 
having the described properties. The air forces used are 
those obtained b y a t w o— d i me n s i o na 1 wing treatment, Good 
judgment • in selecting the representative wing section 
gives calculated critical flutter speeds that come close • 
to e xpe r i ae nt r> 1 v a lue s • 

When one or both of the coupled nodes of vibration 
making up the flutter vary in magnitude along the span or 
do not involve t he entire wing, the effect of t he wing 
shape in flutter is introduced into the problem by apply- 
ing we ight ing factors to t he c oupl ing terms present in 
the equation of notion of the f lut t er ing wing. F or two- 
dimensional -flutter cases of an inf in it e wing involving 
the aileron, the weighting factor t o be applied to the 
a i le r on— bending and ai ler on— t or s i on c oupl ing t firms w ould 
be the ratio of aileron spnn to wing span. For an actual 
wing in the case of b e nd ing-t or s i on flutter, the weighting 
factor £ would depend on the wing bending and wing tor- 
sion shapes in the flutter vibration. 

If th<* displacement of any point of the wing due to 

the bending mode of vibration is given by h (x ) e 1 (^t+cp i ) 
and the displacement due to torsional vibration is written 
£S a(x)e i + *i) 

where 

x span coordinate measured from wing root 
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h(x) 



amplitude of bonding vibration at x 



a (x) 



amplitude of torsional vibration at x 



flutter frequency 



t 



t ime 



Vai^'i phas<* angles of the respective vibrations 

Then, for the ce.se of a rectangular wing with section 
parameters constant along the span, the coefficient of 



whoro I is the wing semispan length* 



THE US£ OP LAGRANGE ! S EQUATIONS IN* TES FLUTTER PROBLEM 



5? ho application of Lagrange's equations and general- 
ized coordinates to the flutter problem provides a means 
of treating the wing as a continuous structure capable of 
vibrating in its various modes with amplitude variable 
along the span* The characteristic determinant developed 
in the course of the solution of a flutter problem yields 
upon the solution of the equations derived therefrom an 
estimate of the shape assumed by the fluttering wing as 
well as the critical flutter speed and vibration frequency. 
One weakness of this method at present is the practical 
necessity for assigning to each wing section the air forces 
that would exist on a two-dimensional wing vibrating in a 
manner identical with the section, This assumption leaves 
out the effect of aspect ratio. For the higher reduced 
flutter frequencies the error introduced by neglecting the 
effect of aspect ratio is not large. 

In applying Lagrange's equations it is considered that 
the function which, gives the flutter vibration shape along 
the span for a given case can be made up of a sum of ele- 
mental span functions, each multiplied by a factor. These 
multiplying factors are taken as the generalized coordi- 
nates of the problem and their evaluation gives the func- 
tion expressing the variation of vibration amplitude along 
the wing span for each vibration mode. If H(x) is the 
amplitude span function of the bending vibration and a(x) 



coupling is obtained by the integral 




o 
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4 

to 

rH 
J 



the corresponding function for. the torsion vibration, 
then 



H(: 



- i (u)t + cp n ) v 
:) = Sq n o TI1 h n (x) 



and 



au) = Ep n e n a n (x) 



g , p n gonoraliz^d coordinates 

h n (x),a n (x) elemental amplitude functions for "bonding 

and torsion vibrations 

9 n , phase angles of banding and torsion vibration 

nodes 

If the wing bending shape in flutter is assumed to 
be some combination of th<* wing shapes in th--> first and 
second natural., vibr-at i on modes- of the wing., t.he use of 
two oleir^-ntal functions h x (x) and h 5 (x) , representing 
the first and. second natural vibrat ion" bend ing ;v.odes, is 
suff ic ie nt f or the s or ies- expre s s ing the flutter binding 
shape- E(x). Because the lowest nature! torsion fre- 
quency of ;the .wing -is so., pny t i:^s. , higher .than the bend- 
ing frequency, the f irs t natural torsion mode of the wing 
is assumed to be a good approximation to a(x), the tor- 
sion flutter sh'.pe, Th* bending and torsion flutter am- 
plitude functions employed are 

\ i(o)t + cp 1 ) l f N " i(wt + cp 2 ) 

Mx) = q x e Yj hi(x) + q g c h 2 (x) 

a(x) = $ x e v Y1/ a 1 (x) 

For t,he case of a rectenguler wing with constant 
section parameters, the various bending find torsion nat- 
ural vibration modes are well known and c- c in be approximated 
ty simple polynomials in x. There is then for h 1 (x), 
hs(x) , and a 1 (x) 
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h a (x) = 12 



1 2 x 2 - - + i x H- 
\ 3 3 J ■ 



h 3 (x) = (1.6276 x 2 - 4.1621 x 3 + 3.3484 x 4 - 0, 9231 3c 5 )l7 2 
ai (x) = (2 x - x 2 ) H 3 

wher" x = — , * being th<" wing seaispan length, and IT 

V 

with subscripts represents normalizing factors so chosen 

pi h 8 pi 
that / ---(x)dx = 1 -/ a 2 (x)dx. Th* values of S 
4 b 3 Jo 

thus d^te rained are 

IT , = 1 . 9?b 
F s = 19.3b 



!I 3 = 1.37 



where b is. the wing seinichord length. Th~ shapes of 
thes« functions ar--> illustrated in figur - ' 1. 

For the b end ing— t or s i on flutter case for the uniform 
rectangular wing, the application of Lagrange's equations 
giv^s for the generalized equations of motion: 

I I _ 

•A / g / \ , 1 OrJ f fb<x x (x)V a . 
! A aa/ aj 3 (x)dx + — - / ( U ' dx pj. 

u J *u> 2 Kb "J \ dx / j 



^ 1 !*A h r^ ~ 

a (x) a a (x)dx I <ii + hg(x) aj(x)dx | q 2 = 0 



*r-hP" 

L "7 h 



J L b J 0 " J 
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H 



r 

| A ca / h x (x)a x (x).dx p x 
u J 0 J 

+ / h, x dx + — - — / I l dx a x 

+ !---/ h 1 (x)h 2 (x)dx+ — / — --i- — d 

i_ b J K*^ KTj/- Ox 2 bx s 



!A C a| h s (x) a x (x)dx ip a 



kau^ M.fc. / ax 2 

0 



i en, g / \ i sir /d s h 3 (x)\ . 
+ !-— h 3 *(x) dx + — - -- / ! : ax 



dx 



<li 



as = 0 



The symbols used ?re consistent with those pcployod 
in r n f rt rencps I -hd 3. For ref ornnc p , the terms ar^ 
1 i s t p d p. s 

M w i ng was s por unit 1 - ngt h 
"fa wingseiaichord 

SI wing bending stiffness; product of Young's modulus 
and E-om^nt of inertia of wing cross section 

GJ wing torsion stiffness 

A pa , etc. complex terns defined in references 1 and 2 

P° rf or i^ing the integrations indicated f or the t e r ias 
containing wing structural stiffness terms yields- 



I 2 

SI f /S 3 h 1 (:c)\ „ ■ . 

— / 1 — t~ — ; *r. - t»h, 1 c 



o 

,1 , 2 



Mo J V ox / h * 



0 



Tho regaining integrals yield 
I 

h z (x) :i a (x) * -0.126 I b £ 

'0 
4 



1 



h A (x) a x (x) dx - 0. 95 6 b 



o 

r 



/ h a (x) a x (;x) dx = 0.172 1 b 



wh^pp 



1 



first natural wing "binding frequency 



^hs second natural wing "bonding fr^quoncy 
first naturrl wing torsion frequency 
critical flutter speed 



uu 



critical flutter frequency 
k poduc- d flutter frequency ~ 

r a radius of gyration of wing mass about static torsion 
axis 

p air density 
and 
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H 



x = 



f VaV /IN 



M 

. The ge n^r ? 1 i ze d equations of motion become, aft or 
s om<- r -duct i on , 

(A aa + X) Pl + (0.95 A ah )q x + (0.172 A ah )q 2 = C 

(0,95 A ca )Pi-+ (A ch f Q 1 X)q 1 + (-0.126 A ch + 0 . 010 95 Q 2 X ) q 2 = 0 
(C.172 A ca ) Pl + (-0.126 A ch + 0.4365^ 1 X)q 1 +^A ch + Q 2 £)q 2 = 0 

The critical values of X and k ar<=» obtained by the 
s i-nultano ous solution of th* r-al ?nd imaginary, equations 
in X and k resulting fro:., the expansion of the chnrac- 
teristic determinant forced from the coefficients of the 
gen^r ali ze-d coordinates of this s-e t of equations-. The 
flutter sp^ed is then obtained from the definitions of these 
tor as 

v / K k v X 

With v and k known, the flutter frequency uu is 
obtained fro:a the expression defining k« 

When the values of the vr-rious critical flutter con- 
stants are substituted in the equations of motion, the 
generalized coordinates are evaluated and the approximation 
to the wing flutter shape is obtained. 
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COMPARATIVE CALCULATIONS 01 CRITICAL FLUTTER SPEEDS 



Critical flutter speeds have "been determined for nine 
rectangular wings in the NACA 8-foot high-speed wind tun- 
nel at LKAL # Por each of these wings six theoretical crit- 
ical speed determinations were made: 

1, Lagrange^ equation method for one "bending and 
one t ors ion mode 

2. Lagrange's equation method for two "bending and 
one torsion mode, with h 1 (x) h 2 (x) dx = -0.125 I = % 

S. Lagrange's equation method for two tending and one 

/ * 

t or s i on mode , with ^ h 1 (x ) h s (:: ) dx arb i trar ily set 
equal to zero (t » 0-) • The coefficient ~0 P 126 I comes 

ri 

from the term J o h x (x) h 2 (x) dx which would he sero if 

h 1 (x) and h 2 (x) were a pair of strictly orthogonal 
f unct i oris • 

4« Lagrange's equation method v/hen ono torsion mode 
and tv. r o "bending modes are used; the second mode is a tip 
deflection made up of that portion of h g (x) beyond the 
node occurring a o out C ,71 along the semispan. These com- 
putations are given for plastic wings, which showed a tip- 
deflection tendency in t he w ind— tunnel tests . 

5. The odor sen 1 s method when the first natural "bending 
frequency is used 

6 # Iheodorsen's method when the second natural "bend- 
ing frequency is used 

The Justification for making calculations with I 
equa 1 to zero (gr oup 5 ) is based on the argument that the 
small change in the shapd of the second "bending mode nec- 
essary to make it orthogonal with the first bending mode 
is so small that the other coupling terms involving the 
second bending mode are not changed appreciably. Previous 
computations have shown that small changes in these coup- 
ling terms do not produce an appreciable effect on the com- 
puted critical flutter speed. 

The critical-flutter-speed computations with 



Theodorsen Vs method using the second "bonding frequency 
were made to show that the experimental critical, flutter 
speed usually lies between the computed values obtained 
with the first and second wing— "bending frequencies. This 
point is discussed in do ta i 1 in reference 2 . 

The values of the wing parameters used in these 
computations are given in table I. The results of these 
computations, as well as the experimental critical flutter 
speeds for the corresponding wings, are given in table II. 

ZProrn a comparison of the values obtained for the 
critical flutter speeds computed by Theodorsen ! s method 
and Lagrange's equation method, it appears that considera- 
tions involving the effect of wing flutter shape in the 
computation of critical flutter speeds may not he neces- 
sary in view of the small and inconsistent differences in 
the v c, 1 u o s o b t a i n e d a n d the increased a m o u n t o f w ork in— 
tr educed, 

tfings 3 to 14A were of largo aspect ratio and tnere- 
fore were suitable for introducing additional bending modes 
into the shape of the fluttering wing. ?or those wings 
Lagrange's equation ..method yields results that do not rep- 
resent a consistent improvement in the predicted critical 
flutter speeds as compared with the values obtained by 
The odor con's method. 

J or w ings 16 to 18, which were of plas tic mater ial 
constx uc ted more like actual wings , t ho value s of the 
critical flutter speeds obtained with Lagrange's equation 
methoa using two bending modes, in general, come further 
from the experimental values than do The odor s en ' s . Per 
these wings, the second bending mode introduced into the 
p r o b 1 e m a ~o p e a r s t o h a v c no e f feet on t h e c a 1 c u i a t e d crit- 
ical flutter speed. Because these wings are weak in bend- 
ing near the tip, choosing the second bending mode as a 
tip deflection tends to bring Lagrange's equation values 
a little closer to the exper imo nt al flutter speeds. The 
t ip-dof loot ion function was taken by using the expr e s s i on 
h g (x) and f ixing its value as zor o up t o the pas ition of 
the node that occurs about 0.71 along the span from the 
root. This effect emphasises the point that in limiting 
La gr a n g e ' s e a u a t i on m e t h o d t o a s ma 1 i nu mb e r of e 1 e m o n t a 1 
shaoos , a great deal depends on the proper choice of these 
funct i ons • 



¥hen one bending and one torsion mode arc used in 
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Lagrange *s equation method, the results are, in general, 
less accurate than when a second heading node is intro— 



The application of Lagrange *s equations to the flut- 
ter problem described in this report is similar to that 
developed in reference 3 by S , J. Lorin^j, In carrying 
out the computations described in this section of the 
rep or t t he rallies sub s t it ut ed for t he fro q v e ncy t er lis in 
some of the flutter calculations were obtained from wing— 
vibrat i on measur e mont s w-hi ch ^vo c ouplo d f r equenc io s 
that may differ from the true mode frequency by. as much 
as 10 percent • In this respect the computations with 
Lagrange's equations n-ado in this report differ somewhat 
from the method outlined by Lor in g in which uncoupled val- 
ues of the wing frequencies arc specified. In comparison 
with the errors introduced by uncertainties regarding the 
location of the static torsion axis and particularly by 
neglecting v ibra t i on damping , the err or in c omput o d crit- 
ical flutter speed resulting from the use of coupled 
vibration frequencies is not largo • 



!• Vrhen limited by practical considerations to the 
use of two bending— mode functions , Lagrange ! s equation method 
gives predict od critical flutter s ponds that are close to 
the values obtained by application of Iheodorsen's flut- 
ter theory* There is no indication from the values ob- 
tained by cither method in this study that one method 
will consistently give predicted flutter speeds closer to 
the oroor imental values than will the other • 

2 # J or wings of normal aspect ratio, the introduction 
cf a second bonding mode in Lagrange T s equal ion system 
appears to have little effect on the' critical flutter 
speeds obtained, Creed judgment must be used in selecting 
the elemental function representing the second bending 
flutter mode if an appreciable effect is to bo made on 
the calculated critical flutter speed with a second bend- 
ing mode introduced into the problem. 

Langley Homer ial Aeronautical Laboratory, 

national Advisory Commit teo for Aeronautics, 
Langley Field, 7a. 
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zero reduces th.o accuracy of the rpsults. 



equal to 



CONCLUSIONS 
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TABLE I.- WING PARAMETERS USED IN FLUTTER COMPUTATIONS 



Wing "a ^ * K a r a 2 

(cpro) ' w a UJa (ft) 



2 


1085 


0.0694 


0,U25 


0.^00 


6 


1150 


.120 


.670 


.5012 


10 


1470 


.0694 


.425 


.5025 


11 


1790 


.0659 


.405 


.500 


12 


1925 


.0675 


.411 


.502 


13 


2350 


.0642 


.400 


.500 


lU 


. 2155 


.0659 


.399 


.5013 


16 


730 


.1192 


.333 


.500 


16s 


795 


.132 


.906 


.500 


17 


1600 


.247 


1.50 


.500 


is 


1650 


.252 


1.612 


.507 



' 1/90 


-0.400 


0.25 


O.3I25 


0.1170 


-.1762 


.03 


.272 


.0113S 


-.400 


.25 


.3oo 


.01052 


-.400 


.254 


.3045 


.00972 


-.251 


.03 


.243 


.00673 


-.too 


.2316 


.2385 


.00772 
.0746 


-.39s 


.234 


.2397 


-.450 
-.U50 


.093 


.36O 


.0743 


.093 


.360 


.0349 


-.104 


.052 


.345 


.0796 


-.450 


.093 


.360 
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